Quantum noise in three-dimensional BEC interferometry 
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We develop a theory of quantum fluctuations and squeezing in a three-dimensional Bose-Einstein 
condensate atom interferometer with nonlinear losses. We use stochastic equations in a truncated 
Wigner representation to treat quantum noise. Our approach includes the multi-mode spatial evolu- 
tion of spinor components and describes the many-body dynamics of a mesoscopic quantum system. 
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Atom interferometry is an important quantum tech- 
nology at the heart of many proposed future applica- 
tions of ultra-cold atomic physics. Bose-Einstein con- 
densates (BECs) or atom lasers are macroscopic quan- 
tum objects and have potential advantages as interfero- 
metric detectors and sensors, provided one can precisely 
extract atomic phase information. However, unlike pho- 
tons, atoms can interact strongly causing dephasing and 
loss of interference fringes. An intimate understanding of 
quantum many-body dynamics is the key to calculating 
interaction-induced dephasing in the measurement pro- 
cess. This is essential for a quantitative theory of atom 
interferometry. 

In this Letter we present a simple, yet quantitatively 
accurate theoretical approach to simulating the dynam- 
ics and evaluating limits of atom interferometry at large 
atom number, using a truncated Wigner representa- 
tion [lHl]. This method extends the conventional Gross- 
Pitaevskii equations describing a Bose condensate to in- 
clude quantum noise effects, including noise due to linear 
and nonlinear losses. The theory allows the accurate in- 
clusion of quantum fluctuations due to nonlinear losses, 
which is a dominant effect when atom numbers are in- 
creased to improve fringe visibility. 

Importantly, we can clearly demonstrate where fringe 
visibility is driven by quantum fluctuations, and where it 
is driven by trap inhomogeneity and dynamical effects, 
in order to choose optimal conditions for quantum noise 
reduction and spin squeezing. These calculations are a 
first step towards understanding mesoscopic superposi- 
tions and entanglement in ultra-cold atomic gases. An 
advantage of our method compared to the variational ap- 
proaches used elsewhere 0, HI is that it allows us to treat 
a large number of independent field modes and particles, 
thus including degrees of freedom that are excited due 
to collisional and nonlinear loss dynamics d 0. Our 
theory can be readily extended to include finite temper- 
ature initial conditions .2, 8], which will be treated else- 
where. Nonlinear losses and finite temperature effects 
can be also described within the confines of the varia- 
tional approach ■ 

Quantum phase-diffusion is defined as the phase noise 



induced by number fluctuations which are conjugate to 
phase. This is a fundamental feature of BEC interferom- 
etry, and can only be removed when there are no interac- 
tions. However, there are other reasons for decoherence, 
which are also important. The approach used here cap- 
tures all three significant features of atom interferometry 
that can result in decoherence: phase-diffusion, losses, 
and trap inhomogeneity effects. The results given in 
this paper are applicable to simulations where the atom 
number per lattice point or mode is large. We focus on 
trapped BEC interferometry, for definiteness, although 
the method has a more general applicability. 

We start by assuming that the BEC has s-wave interac- 
tions, together with Markovian losses due to n-body col- 
lisions. We employ a master equation together with the 
Wigner-Moyal quantum phase-space representation [llj 
and a truncation of third and higher-order derivatives in 
the equations of motion. If we regard the commonly used 
Gross-Pitaevskii equation as a classical, first approxima- 
tion to mean- field condensate dynamics, the truncated 
Wigner approach is best thought of as the second term in 
an expansion in inverse particle number. This truncation 
has been shown to be valid in the limit of large particle 
number [H-Q. It has been further tested by compari- 
son with the exact positive-P simulation method. The 
truncated Wigner method is particularly useful in low- 
dimensional and trap environments, where it has success- 
fully predicted quantum squeezing and phase-diffusion 
effects, in good agreement with dynamical experiments 
in photonic quantum soliton propagation 
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In the present Letter, we treat an ultra-cold, inter- 
acting multi-component spinor Bose gas in D effective 
dimensions. The basic Hamiltonian is easily expressed 
using quantum fields *&j(x) and $!j(x), where *&j(x) cre- 
ates a bosonic atom of spin j at location x, and *&j(x) 
destroys one; the commutators are [^-(x), S^Kcc')] = 
5^ D \x — x')8jk ■ The resulting physics of a dilute, low- 
temperature Bose gas is well-described in the s-wave scat- 
tering limit by an effective Hamiltonian with contact in- 
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teractions and external potentials: 



U 



H/H = I d u x I *}K jk * k + ^$t$t^ fe ^ 



(1) 



Here we omit the field argument (x) for brevity, and 
use the Einstein summation convention of summing over 
repeated indices. Kjk is the single-particle Hamiltonian: 



Kjk = r 2^ v + Ui + v ^/ h ) 6 > k + n »®> (2) 

where m is the atomic mass, V} is the external trapping 
potential for spin j, ujj is the internal energy of spin j, 
fijfc represents a time-dependent coupling that is used to 
rotate one spin projection into another, and Ujk is the 



atom-atom interaction term. Thus, 



= is the 



3 ~ 

spin-j atomic density. For a dilute gas at low enough 
temperatures, Ujk = -iirhajk/m, where ajk is the s-wave 
scattering length in three dimensions. Here we assume 
a momentum cutoff k c <C 1/ajk, otherwise the couplings 
must be renormalized [3j. 

We proceed by using a stochastic phase-space method 
that allows a numerical simulation of the quantum dy- 
namics 0, 0, EH- Defining a Wigner function W(^f), 
where VP is a c-number field corresponding to the quan- 
tum field "P, this has a unitary time-evolution equation: 



dW 
~dt 
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Next, higher derivative terms of type O [5 3 /(5\E , |] are 
truncated. This approximation neglects higher-order 
terms in an expansion in 1/y/W, and is therefore valid 
in the limit of N 3> M where N is the atom number and 
M is the number of low-energy modes included [H, 03, [H| . 
In free-space calculations it is important to maintain this 
mode truncation. In the relevant limits where the tech- 
nique is applicable, the equations simply reduce to Gross- 
Pitaevskii equations with Gaussian fluctuations of the 
initial conditions: 



—± = -i {Kjk^k + U jk \^k 12 



(4) 



For initial conditions in interferometry it is usually suf- 
ficient to consider a coherent state amplitude corre- 
sponding to a typical initial state with Poissonian num- 
ber fluctuations, as produced by a beam-splitter. In this 
case the initial Wigner amplitude has a Gaussian random 
distribution, with ^j(x,t ) = *f(sc) + A$Fj(x), where: 
(A^ j (x)A^* k (x')) = S jk S D (x - x')/2. This initial noise 
is necessary because the Wigner representation generates 
symmetrically ordered correlation functions, and includes 
vacuum fluctuations. For greater accuracy, the initial 
state can be modified to account for initial correlations, 
thermal noise, or additional fluctuations. If normal or- 
dered correlations are measured, one has to express them 
as a sum of symmetrically ordered terms. 



This includes all the known nonlinear quantum noise 
effects of quantum dynamics, like phase diffusion, entan- 
glement and quantum squeezing, in the limit of large par- 
ticle number. The initial noise terms do not occur in the 
semi-classical Gross-Pitaevskii approximation, which is 
therefore unable to predict these effects. Thus, while the 
lossless equations are identical to the Gross-Pitaevskii 
equations, the inclusion of initial noise terms together 
with nonlinear interactions leads to quantum phase- 
diffusion. Such methods can be used for either free-space 
or trapped atom interferometry, provided there is an ap- 
propriate mode truncation. 

Additional quantum noise enters from the effects of 
damping and losses, due to the fluctuation-dissipation 
theorem. These effects are important at high densities 
in atomic traps. They can be included via an additional 
Markovian master equation [16| defined so that, 



dp 



H,p 
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d u x4 n) [ P ] , (5) 



where n is the number of interacting particles, I = 
(7i, I2, ■ ■ ■ , In) is a vector indicating the spins that are cou- 
pled, and we have introduced local Liouville loss terms, 



4 n) [p] = 2dl n) pO^ 



d^d^p- po\' l,l o\ n> . (6) 
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The reservoir coupling operators O t are the distinct n- 

2 (n) 

fold products of local field annihilation operators, 0\ — 
6[ n \$) = ^!i 1 (x)^i 2 (x) . . . $;„(a;), describing local n- 
body collision losses. 

After transforming these new terms to evolution equa- 
tions for the Wigner distribution, the drift term Aj 
changes the Gross-Pitaevskii evolution to include non- 
linear damping, while the next terms in the evolution 
equation give rise to additional Fokker-Planck diffusion 
terms associated with quantum noise from the loss reser- 
voirs, given by: 



j x _ An) aol n) * doj n) 

dm d^ k 



W. (7) 



This leads to a stochastic equation: 



dt 



n . I 

(8) 

where the nonlinear loss has the form: 



111 



and { x ->t) is a corresponding complex, stochastic 
delta-correlated Gaussian noise with 

(C[ n \x, t)ci m) *(x>, f )) = 8 lh 5 nm 5 D (x x') 5 (t ~ f ) . 

(10) 



3 



The multiplicative noise coefficient 



(n) 



(11) 



is a fluctuation-dissipation term, so that the Wigner vari- 
ables remain equivalent to the corresponding operators. 

The loss coefficients in equations ([5]) can be converted 
to the conventional form, which is defined using atom 
number losses: 



(n) mi rri2 
n 3 = -%j n l n 2 



(12) 



where rij is the density of component j and rrij is the 
number of spin-j atoms lost in the collision. The conver- 
sion can be carried out as — 2m jK^ 1 . 

In this work we use a basis of plane waves in the volume 
V, and the density of component j is calculated as a 
probabilistic average: 



nj(x) = (**(x)* j (a;))p ath8 



M 

w 



(13) 



Here we use the fact that the approximate Wigner func- 
tion is a probability distribution equivalent to an aver- 
aged sum over different simulation paths. 

To illustrate the applications of this method we 
consider recent interferometry experiments with a 
two-component BEC involving two hyperfine states 
\F = l,m F = -1) and \F = 2,m F = +1) in 87 Rb [l7|. 
A conventional Ramsey sequence (Fig. [TJ (a)) has been 
used with a BEC confined in a cigar-shaped magnetic 
trap with the frequencies (97.0, 97.0, 11.69) Hz in a bias 
magnetic field of 3.23 G, so that magnetic field dephas- 
ing is largely eliminated [l3|. The first ir/2 pulse pre- 
pares a non-equilibrium superposition of states |1,— 1) 
and |2, +1) and the spatial modes of two components pe- 
riodically separate and merge again [l8|- The spatially- 
separated spin components evolve differently, as they 
have different scattering lengths. As a result, these col- 
lective oscillations lead to periodic dephasing and self- 
rephasing of the BEC components, clearly visible in both 
GPE and Wigner simulations of interference fringe vis- 
ibility V (Fig. [TJ (b)). Asymmetric losses of two states 
are one cause of the contrast decay. This can be par- 
tially compensated by the application of a spin echo pulse 
mid- way through the evolution (Fig. (TJ (c)). The GPE 
simulations wrongly predict (dashed lines) that visibil- 
ity is largely recovered at long evolution times using the 
spin echo method. However, the addition of quantum 
noise (solid line) via the Wigner simulations noticeably 
speeds up the visibility decay even with a spin echo pulse 
present. This is in agreement with experimental observa- 
tions, and shows that these effects play a significant part 
in the decay of visibility, even for large particle numbers. 

The important feature of these quantum dynamical 
simulations is that they are able to treat large numbers 
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FIG. 1. (color online) Timeline of the experiment for Ram- 
sey (a) and Ramsey with spin echo (c), (b) and (d) are the 
simulated plots of interferometric visibility. Classical GPE 
(red dashed lines) and Wigner calculations (blue solid lines) 
are shown. N = 5.5 x 10 4 , u x = w a = 2tt x 97.0Hz, 
cj z = 2tt x 11.69Hz, an = 100.4 a , a 12 = 97.993 a , a 22 = 
95.57 ao [17| |. ao is the Bohr radius. Nonlinear atomic losses: 



Tifi = 5.4 x 
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= 1.51 x 10" 



of atoms (55,000 in this case), while correctly tracking 
all the quantum noise sources, and also extending the 
simulations to long time-scales. Both of these features, 
large atom numbers and long time-scales, are essential in- 
gredients to accurate interferometric measurements. The 
simulations give accurate predictions despite large, multi- 
mode dynamical motion in three dimensions and sub- 
stantial losses of most of the condensate atoms Uj. On 
longer time-scales, the experimental accuracy is limited 
by technical noises, and we have no data for comparisons. 

The Wigner method is able to predict not only the 
average values of different observables, but their vari- 
ances too. As an example, we now calculate the de- 
gree of squeezing £ 2 , introduced in 2l|, l22| . This allows 
the analysis of experiments beyond the usual shot noise 
limit 23j, |24| . If we define spin component operators as 
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FIG. 2. (color online) Wigner simulations of squeezing in 
the vicinity of 9.1 G Feshbach resonance in 87 Rb. Simulation 
parameters are the same as in Fig. [TJ except for the inter- 
component scattering length and corresponding loss rate [2(j : 
ai2 = 80.0 a , j[f = 3.88 x 10~ 12 cm 3 /s (blue solid line), 
ai2 = 85.0 a , 7^ = 1-95 x 10~ 12 cm 3 /s (red dashed line), 
ai2 = 90.0 do, 712' = 7.13 x 10~ 13 cm 3 /s (green dash-dotted 
line) and a 12 = 95.0 a , y$ = 8.54 x 10" 14 cm 3 /s (black 
dotted line) 



S x = [S+ + S\}/2, S y = i[S+ - SU/2 and 



5+ = / d u x 



*t(x)*i(»)-*];(aj)*2(a!) 



(14) 



the average values of spin vector components can be cal- 
culated similarly to the average density in equation (|13[) . 
Now we can shift the centre of coordinates to the end of 
the vector (S) = {(S x ), (S y ), (S z )} and rotate it, making 
the z' axis parallel to S and choosing the remaining axes 
x' and y' so that ASy is minimized. 

Second-order moments of spin operators are repre- 
sented using fourth-order moments of field operators 
and, therefore, can be calculated using wavefunctions in 
Wigner representation. In this new coordinate system 
the squeezing parameter is expressed simply as 



■2 NA ^ 
(5)2 ■ 



(15) 



where N is the total number of atoms. The formula for 
the squeezing parameter in the original coordinate sys- 
tem can be found elsewhere [i[ . The squeezing parameter 
serves as a measure of atomic entanglement in the con- 
densate, when the atoms are entangled if £ 2 < 1 0. 

As an illustration of the power of the Wigner 
method and of the effect of two-body losses we con- 
sider the temporal evolution of the squeezing parameter 
in a Ramsey interferometer for optically trapped two- 
component BEC in 87 Rb (states \F = 1, rap = +1} and 
\F = 2, rap — —1)) near a Feshbach resonance at 9.1 G 
(Fig. [2]). Different inter-component scattering lengths 
(corresponding to different magnetic fields) were used in 
order to find the optimal regime with significant squeez- 
ing. The value a%2 — 80ao provides stronger nonlinear 



interactions, but larger two-body losses quickly eliminate 
the squeezing effect. Feshbach tuning to a\i = 90ao en- 
sures the best squeezing (— 6dB at 40 ms), whereas long 
lasting squeezing is predicted for a%2 — 95ao- 

These simulations predict the degree of quantum noise- 
reduction, which is a technique that can be used to im- 
prove precision quantum interferometry. The quantum 
noise is initially reduced due to a stretching and rota- 
tion of the quantum noise elli pse , similar to that found 
in quantum soliton squeezing 12, 25|. For lossless envi- 
ronments, this effect is optimized when the cross-species 
scattering length is maximally different to the intra- 
species scattering length, which is found near the Fesh- 
bach resonance. However, there are competing nonlinear 
loss effects at the Feshbach resonance, which means that 
some detuning of the magnetic field is essential to re- 
duce these detrimental losses. Subsequently, the squeez- 
ing is destroyed in time as the two fields recombine and 
interfere with each other. Importantly, these quantum 
squeezing calculations indicate conditions that will allow 
this macroscopic quantum effect to be experimentally ob- 
served in ultra-cold atomic BEC for much larger atom 
numbers than calculated previously. 

In conclusion, the truncated Wigner- Moyal quantum 
phase-space theory is well-suited for a variety of cal- 
culations of many-body dynamics of quantum degener- 
ate gases. Although the Wigner function is not a true 
probability, it becomes one to a good approximation 
in the limit of large particle number, so that a prob- 
abilistic (random) sampling of phase-space trajectories 
can be used to treat quantum fluctuations. Further- 
more, it is straightforward to treat a spatially inhomo- 
geneous, many-mode problem as found in an extended 
3-dimcnsional environment. Initial conditions such as a 
coherent beam-splitter as employed in typical atom in- 
terferometry experiments can be readily included, with 
proper accounting for both quantum effects and nonlin- 
ear losses in the condensate. 
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